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o ' Abstract: The production of four bottom quarks is an important benchmark channel 

^ I for Higgs analyses and searches for new physics at the LHC. We report on the calculation 

of the next-to-leading order QCD corrections to the process pp — )• hhhh + X with the 

Helac-NLO automated framework, and present results for inclusive cross sections and 

K* ■ differential distributions. We discuss the impact of the higher-order corrections and, in 

vQ . particular, the effect of the bottom quark mass. In addition, we provide an estimate of the 

OO , theoretical uncertainty from the variation of the renormalisation and factorisation scales 

and the parton distribution functions. The results are obtained with a new subtraction 

formalism for real radiation at next-to-leading order, implemented in the Helac-Dipoles 
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1 Introduction 

The production of four bottom quarks, pp — >• bbbb + X, is an important background to 
various Higgs analyses and new physics searches at the LHC, including for example Higgs- 
boson pair production in two-Higgs doublet models at large tan/3 [1], or so-called hidden 
valley scenarios where additional gauge bosons can decay into bottom quarks [2]. Accurate 
theoretical predictions for the Standard Model production of multiple bottom quarks are 
thus mandatory to exploit the potential of the LHC for new physics searches. Furthermore, 
the calculation of the next-to-leading order (NLO) QCD corrections to pp — )• bbbb + X pro- 
vides a substantial technical challenge and requires the development of efficient techniques, 
with a high degree of automation. We have performed a NLO calculation of bbbb production 
at the LHC with the Helac-NLO system [3]. In particular, we present results based on 
a new subtraction formalism [4, 5] for treating real radiation corrections, as implemented 
in the Helac-Dipoles package [6]. Two calculation schemes have been employed, the so- 
called four-flavour scheme (4FS) with only gluons and light-flavour quarks in the proton, 
where massive bottom quarks are produced from gluon splitting at short distances, and the 
five-flavour-scheme (5FS) [7] with massless bottom quarks as partons in the proton. At all 
orders in perturbation theory, the four- and five-flavour schemes are identical, but the way 
of ordering the perturbative expansion is different, and at any finite order the results do 
not match. Comparing the predictions of the two schemes at NLO thus provides a way to 
assess the theoretical uncertainty from unknown higher-order corrections, and to study the 
effect of the bottom mass on the inclusive cross section and on differential distributions. 
First NLO results for pp — t- bbbb + X m. the 5FS have been presented in Ref. [8]. We not 
only provide an independent calculation of this challenging process with a different set of 
methods and tools, but also a systematic study of the bottom quark mass effects by com- 
paring the 5FS and 4FS results. We note that NLO results for the production of four top 
quarks in hadron collisions have been discussed in Ref. [9]. In addition, NLO calculations 



for processes of similar complexity have recently been presented in the literature, including 
NLO QCD corrections to ti production in association with two jets [10], the production of 
a single gauge boson plus jets [11], double gauge boson production with two jets [12] and 
multi-jet production [13]. 

The paper is organised as follows. In section 2 we briefly summarise the set-up of the 
calculation. Numerical results for inclusive cross sections and differential distributions for 
bbbb production at the LHC are presented in section 3. We summarise in section 4. 

2 Theoretical Framework 

The calculation of the process pp — t- bbbb + X at NLO QCD comprises the parton processes 
gg —7- bbbb and qq — )• bbbb at tree-level and including one-loop corrections, as well as the 
tree-level parton processes gg — > bbbb + g, qq ^- bbbb + g, gq ^ bbbb + q and gq — )• bbbb + q. 
In the four-flavour scheme q G {u,d,c,s}, and the bottom quark is treated massive. The 
bottom mass effects are in general suppressed by powers of mf,/fi, where fj, is the hard 
scale of the process, e.g. the transverse momentum of a bottom-jet. Potentially large 
logarithmic corrections oc ln{mf,/fi) could arise from nearly collinear splitting of initial- 
state gluons into bottom quarks, g — )• bb, where the bottom mass acts as a regulator 
of the collinear singularity. This class of ln(7Tife/^)-terms can be summed to all orders in 
perturbation theory by introducing bottom parton densities in the five-flavour scheme. The 
5FS is based on the approximation that the bottom quarks from the gluon splitting are 
produced at small transverse momentum. However, in our calculation we require that all 
four bottom quarks can be experimentally detected, and we thus impose a lower cut on the 
bottom transverse momentum, px^ > Pr™- As a result, up to NLO accuracy the potentially 
large logarithms in the process pp — >• bbbb+X are replaced by ln{mi,/fi) — >• ln(p™™/^), with 
7716 ^ PtT ~ /^' ^^^ *^^ thus much less significant numerically. Therefore, for the process 
at hand, the differences between the 4FS and 5FS calculations with massive and massless 
bottom quarks, respectively, should be moderate, but may not be completely negligible. 

Our calculation is performed with the automated Helac-NLO framework [3], which 
includes Helac-Iloop [14] for the evaluation of the numerators of the loop integrals and 
the rational terms, CutTools [15], which implements the OPP reduction method [16-19] 
to compute one-loop amplitudes, and OneLoop [20] for the evaluation of the scalar inte- 
grals. The singularities for soft and collinear parton emission are treated using subtraction 
schemes as implemented in Helac-Dipoles [6], see the discussion below. The phase space 
integration is performed with the help of the Monte Carlo generators Helac-Phegas [21- 
23] and Kaleu [24], including Parni [25] for the importance sampling. 

The Helac-Dipoles package has been based on the standard Catani-Seymour dipole 
subtraction formalism [26, 27]. We have now extended Helac-Dipoles by implementing 
a new subtraction scheme [4, 5] using the momentum mapping and the splitting functions 
derived in the context of an improved parton shower formulation by Nagy and Soper [28] . 
Compared to standard dipole subtraction, the new scheme features a significantly smaller 
number of subtraction terms and facilitates the matching of NLO calculations with parton 
showers including quantum interference. The results presented here constitute the first 



application of the Nagy-Soper subtraction scheme for a 2 — >• 4 scattering process with 
massive and massless fermions. A detailed description of the implementation of the new 
scheme, and a comparative study of the numerical efficiency and the speed will be presented 
elsewhere. 

3 Numerical Results for the LHC 

In this section we present cross-section predictions for the process pp — >■ bbbb + X at the 
LHC at the centre-of-mass energy of y/s = 14 TeV. We discuss the impact of the NLO-QCD 
corrections, and study the dependence of the results on the bottom quark mass. We show 
results obtained with the standard dipole subtraction scheme for real radiation, and with 
the new subtraction formalism [4, 5] as implemented in Helac-Dipoles. 

Let us first specify the input parameters and the acceptance cuts we impose. The 
top quark mass, which appears in the loop corrections, is set to mt = 173.5 GeV [29]. 
We combine collinear final-state partons with pseudo-rapidity |ry| < 5 into jets according 
to the anti-fc^^ algorithm [30] with separation R = 0.4. The bottom-jets have to pass 
the transverse momentum and rapidity cuts pT,b > 30 GeV and \yh\ < 2.5, respectively. 
The renormalisation and factorisation scales are set to the scalar sum of the bottom- 
jet transverse masses, fiR = fip = fJ-o = Ht, with Ht = friTfi + fnTl + ''^T.h + "^t5 
and the transverse mass rriT.b = \ h^t +Prfe- ^^"^ ^^^ five-flavour scheme calculation with 
massless bottom quarks the transverse mass equals the transverse momentum, m^^ b = pT,b- 
Note that the implementation of a dynamical scale requires a certain amount of care, as 
the subtraction terms for real radiation have to be evaluated with a different kinematical 
configuration specified by the momentum mapping of the subtraction scheme. Comparing 
our results as obtained with the Catani-Seymour subtraction and the Nagy-Soper scheme, 
which is based on a different momentum mapping, provides an important and highly non- 
trivial internal check of our calculation. 

3.1 Massless bottom quarks within the five-flavour scheme 

Results are presented for the NLO CTIO [31] and MSTW2008 [32] parton distribution 
functions (pdfs) with five active fiavours and the corresponding two- loop Og. To study the 
impact of the higher-order corrections, we also show leading-order results obtained using 
the CT09MC1 [33] and MSTW2008 LO pdf sets and one-loop running for a^. 

We first discuss the impact of the bottom-quark induced processes, bb — t- bbbb, on the 
hadronic cross section at leading-order. The difference between the qq initiated processes, 
qq — )• bbbb, with and without bottom-quarks is at the level of 2.5%. Moreover, at the central 
scale, (Ji, = Ht, the hadronic cross section is completely dominated by gluon-fusion, with 
only about 1% contribution of all quark-antiquark annihilation processes. The bottom- 
induced contributions are thus negligible, and we decided to neglect bottom initial states in 
the computation of the cross section both at LO and NLO. Note that the suppression of the 
bottom-induced processes, which include for example potentially large forward scattering 
of bottom-quarks through i-channel gluon exchange, depends crucially on the transverse 
momentum cuts we impose on the bottom-jets. 
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Figure 1. Scale dependence of the 5FS LO and NLO cross sections for pp — > bbbb + X at the 
LHC (\/s = 14 TeV). The renormalisation and factorisation scales are set to a common value 
fifi = fip ~ S, fJ-o where /Xq ~ Ht- The CT09MC1 and CTIO pdf sets have been used for the LO and 
NLO cross sections, respectively. 



Let us first present our results for tlie inclusive cross section pp — )• bbbb + X at the 
LHC (y/s = 14 TeV), including the transverse momentum and rapidity cuts specified 
at the beginning of the section. The NLO QCD corrections strongly reduce the scale 
dependence, as demonstrated in Figure 1. The central cross section predictions are collected 
in Table 1. Varying the renormalisation and factorisation scales simultaneously about the 
central scale by a factor of two, we find a residual scale uncertainty of approximately 30% 
at NLO, a reduction by about a factor of two compared to LO. The size of the -ff- factor, 
K = (Tnlo/clo! strongly depends on the pdf set, with K = 1.15 for CTIO and K = 1.37 
for MSTW2008. We emphasise, however, that the K-factor is an unphysical quantity and 
strongly sensitive to the choice of scale through the large LO scale dependence. 

We observe a difference of about -7% and +11% between the CTIO and MSTW2008 
pdf parametrisations at LO and NLO, respectively. The pdf uncertainty as estimated from 
the MSTW2008 error pdf sets [32] amounts to +7.3% and -1.5% at 68% C.L., and is 
significantly smaller than the scale uncertainty. A more systematic discussion of pdf and 
as uncertainties will thus be referred to a forthcoming publication. 



pp -^ bbbb + X 


c^LO [pb] 


0-NLO [pb] 


-f^ = ctnlo/c^lo 


CT09MC1/CT10 
MSTW2008LO/NLO 


infiq+6i-5(57%) 
J^UO-y_36.4(34%) 

qqq+58.7(59%) 
-^^■^-34.9(35%) 


,„o (,+35.6(29%) 
i/J.D_2gg(22%) 

, oc 7+38.8(28%) 
^"^"•'-30.9(23%) 


1.15 

1.37 



Table 1. 5FS LO and NLO cross sections for pp -^ bbbb + X at the LHC (y/1 = U TeV). 
The renormalisation and factorisation scales have been set to the central value /io = Ht, and the 
uncertainty is estimated by varying both scales simultaneously by a factor two about the central 
scale. Results are shown for the CT09MC1/CT10 and MSTW2008LO/NLO pdf sets. 



An important input for the experimental analyses and the interpretation of the exper- 
imental data are accurate predictions of differential distributions. Our calculation is set 
up as a parton-level Monte Carlo program and thus allows us to predict any infrared-safe 
observable at NLO. Figure 2 shows LO and NLO predictions for the invariant mass of 
the bbbb system (upper left panel), the total transverse energy Ht (upper right panel), 
the transverse momentum of the hardest bottom jet (lower left panel) and the transverse 
momentum of the second hardest bottom jet (lower right panel). We also show the theo- 
retical uncertainty through scale variation and the i('-factor as a function of the kinematic 
variable. It is evident from Figure 2 that the NLO corrections significantly reduce the the- 
oretical uncertainty of the differential distributions, and that the size of the higher-order 
effects depends on the kinematics. For an accurate description of exclusive observables and 
differential distributions it is thus not sufficient to rescale a LO prediction with an inclusive 
ET- factor. 

As discussed in section 2 we have performed the calculation with two different subtrac- 
tion schemes, the standard Catani-Seymour (CS) dipole subtraction, and a new scheme 
based on the splitting functions and momentum mapping of an improved parton shower by 
Nagy and Soper (NS). The comparison between the two schemes for the inclusive 5FS cross 
section is presented in Table 2. For the Catani-Seymour scheme we show results without 
(omax = 1) and with (omax = 0.01) a restriction on the phase space of the subtraction as 
proposed in Ref. [34, 35]. As evident from Table 2, the cross sections obtained using the 
Catani-Seymour (CS) and Nagy-Soper (NS) subtraction schemes agree within the numeri- 
cal uncertainty of the Monte Carlo integration. This result not only provides a validation of 
our implementation of the new subtraction scheme into Helac-Dipoles, but also provides 
a non-trivial internal cross check of the calculation. 

We have also compared the results obtained in the CS and NS subtraction schemes for 
various differential distributions. Some examples are collected in Figure 3. We observe full 
agreement between the predictions calculated with the two schemes. 

3.2 Massive bottom quarks v^ithin the four-flavour scheme 

Within the four-flavour scheme bottom quarks are treated massive and are not included in 
the parton distribution functions of the proton. We define the bottom quark mass in the 
on-shell scheme and use m^ = 4.75 GeV, consistent with the choice made in the MSTW2008 
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Figure 2. Differential cross section for pp -^ bbbb + X at the LHC (y/s = 14 TeV) in the 5FS as 
a function of the invariant mass of the bbbb system (upper left panel), the total transverse energy of 
the system (upper right panel), the transverse momentum of the hardest bottom jet (lower left panel) 
and the transverse momentum of the second hardest bottom jet (lower right panel). The dash-dotted 
(blue) curve corresponds to the LO and the solid (red) curve to the NLO result. The scale choice 
is fj,R = pip = Mo = Ht. The hashed area represents the scale uncertainty, and the lower panels 
display the differential K factor. The cross sections are evaluated with the MSTW2008 pdf sets. 



four-flavour pdf. Tlie central cross section prediction in LO and NLO for ^ = Ht using the 
4FS MSTW2008 [36] pdf are shown in Table 3. Comparing with the 5FS results presented 
in Table 1, we observe that the bottom mass effects decrease the cross section prediction 
by 18% at LO and 16% at NLO. The residual scale dependence at NLO is approximately 
30%, similar to the 5FS calculation. 

The difference between the massless 5FS and the massive 4FS calculations has two 
origins. First, there are genuine bottom mass effects, the size of which depends sensitively 



pp -^ bbbb + X 



csK 
'nlo 



=1) 



[pb] 



a 



CS(a„ 
NLO 



=0.01) 



[pb] 



a 



NS 
NLO 



[pb] 



CTIO 
MSTW2008NLO 



123.6 ± 0.4 

136.7 ±0.3 



124.9 ±0.9 
136.1 ±0.5 



124.8 ±0.3 
137.6 ±0.5 



Table 2. 5^5" NLO cross sections for pp -^ bbbb + X at the LHC (y/s = U TeV). Results are 
shown for two different subtraction schemes, the Catani- Seymour (CS) dipole subtraction, without 
(ttmax = Ij o,nd with (aniax = O.OlJ o restriction on the phase space of the subtraction, and the new 
Nagy-Soper (NS) scheme, including the numerical error from the Monte Carlo integration. The 
renormalisation and factorisation scales have been set to the central value fiQ =^ Ht, and the CTIO 
and MSTW2008NLO pdf sets have been employed. 



pp -> bbbb + X 


o-LO [pb] 


0-NLO [pb] 


K = CTNLO/O'LO 


MSTW2008LO/NLO (4FS) 


^. ^+49.7(59%) 
°^-^-29.6(35%) 


, -I jj „+33.3(28%) 
iJ^o-3_29.0(24%) 


1.40 



Table 3. 4FS LO and NLO cross sections for pp -^ bbbb + X at the LHC (^/^ = U TeV). 
The renormalisation and factorisation scales have been set to the central value jjlq = Ht, and the 
uncertainty is estimated by varying both scales simultaneously by a factor two about the central 
scale. Results are shown for the 4FS MSTW2008LO/NLO pdf sets. 



on the transverse momentum cut. For p™™ = 30 GeV we find a 10% difference between the 
5FS and 4FS from non-singular bottom-mass dependent terms. This difference decreases 
to about 1% for p^™ = 100 GeV. Second, the two calculations involve different pdf sets 
and different corresponding a^. While a 4FS pdf has, in general, a larger gluon flux than a 
5FS pdf, as there is no (7 —t- bb splitting, the corresponding four-flavour as is smaller than 
for flve active flavours. For pp — ?• bbbb ± X the difference in a^ is prevailing and results in 
a further reduction of the 4FS cross section prediction by about 5%. This latter difference 
should be viewed as a scheme dependence rather than a bottom mass effect. 

In Figure 4 we present the differential distribution in the transverse momentum of the 
hardest bottom jet, as calculated in the 5FS with massless bottom quarks and in the 4FS 
with nih = 4.75 GeV. We show the absolute prediction at LO and NLO, and the predictions 
normalised to the corresponding inclusive cross section. The latter plots reveal that the 
difference in the shape of the distributions in the 5FS and the 4FS is very small. We find 
similar results for other differential distributions. 

Let us finally present the comparison of the massive bottom quark results as obtained 
with the Catani-Seymour and Nagy-Soper subtraction schemes, see Table 4. We observe 
full agreement between the two calculations within the numerical error of the Monte Carlo 
integration, and thereby validate our implementation of the NS subtraction scheme also 
for the case of massive fermions. 
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Figure 3. Differential cross section for pp -^ bbbb + X at the LHC (y/s = 14 TeV) in the 5FS as 
a function of the average transverse momentum of the bottom jets (upper left panel), the average 
rapidity of the bottom jets (upper right panel), the total transverse energy (lower left panel) and 
the bbbb invariant mass (lower right panel). The dash-dotted (brown) curve corresponds to the 
Catani- Seymour (CS) and the solid (orange) curve to the Nagy-Soper (NS) subtraction schemes, 
respectively. The lower panels show the ratio of the results within the two schemes. The scale choice 
is fj,fi = lip = Mo = Ht, the cross sections are evaluated with the CTIO pdf set. 



3.3 Comparison witli results presented in the literature 



A detailed comparison of our results with Ref. [8] has been performed. We find agreement 
for the virtual amplitude at one specific phase space point, but cannot reproduce the 
published results for the integrated hadronic LO and NLO cross sections with the setup as 
described in Ref. [8]. When the CTEQ6.5 pdf set [37] is used rather than CTEQ6M [38] 
as specified in [8], and the factorization and renormalization scales are set to the common 
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Figure 4. Differential cross section for pp -^ bbbb + X at the LHC (y/s = H TeV) in the 4FS 
and 5FS as a function of the transverse momentum of the hardest bottom jet. Also shown are the 
normalised distributions at LO (lower left panel) and at NLO (lower right panel). The lower panels 
show the ratio of the results within the two schemes. The scale choice is pL^ = fip = /io = Ht, the 
cross sections are evaluated with the 5FS and ^FS MSTW2008 pdf sets, respectively. 



pp -> bhhh + X 


CS{Omax = l) r 11 

^NLO [pb] 


CS (amax=0.01) r , i 


<lo [pb] 


MSTW2008NLO (4FS) 


118.3 ±0.5 


118.2 ±0.7 


118.0± 0.5 



Table 4. J^FS NLO cross sections for pp -^ bbbb + X at the LHC (y/s = I4 TeV). Results are 
shown for two different subtraction schemes, the Catani- Seymour (CS) dipole subtraction, without 
(ttmax = Ij o,nd with (Ckmax ~ O.Olj restriction on the phase space of the subtraction, and the 
new Nagy-Soper (NS) scheme, including the numerical error from the Monte Carlo integration. 
The renormalisation and factorisation scales have been set to the central value fj,o = Ht, arid the 
MSTW2008NLO 4FS pdf set has been employed. 



value 



f^o 



A\ Z^PT,i 



(3.1) 



the LO result published in [8] can be reproduced as shown in Table 5: 



a 



[8] 
LO 



[pb] 



o-LO [pb] 



94.88 ± 0.14 



94.74 ± 0.20 



Table 5. LO cross section for pp — > bbbb + X at the LHC f^/s — 1^ TeV) in comparison with the 
result of Ref [8]. Results are shown including the numerical error from the Monte Carlo integration. 
The scale choice is pLu = pip = /iq, with pL^ as defined in Eq. (3.1), and the cross sections are 
evaluated with the CTEQ6.5 pdf set. 

After comprehensive numerical checks and the communication with the authors of 

Ref. [8] it turned out that the NLO numbers published in [8] are based on a scale setting 

that mixes partons and jets and that is not consistent with what is specified in Ref. [8]. 

Adopting fiji = fip = hq, with ij,q defined in Eq. (3.1), and summing over the transverse 

momenta of all jets, we obtain the NLO numbers presented in Table 6. Note that at NLO 

the final state can consist of four or five jets, as determined by the jet algorithm. Our 

[81 
results are compared to the NLO number as published in Ref. [8], ctj^Jlq; and to a corrected 

number obtained by means of private communication from the authors of Ref. [8] , Cj^lo 
The corrected result agrees with our calculation. 



f^NLO [pb] 


[8], corr. r ,-, 
^NLO [pb] 


CS (amax=0.01) r ni 

f^NLO [pb] 


^NLO [pb] 


<0 [pb] 


140.48 ± 0.64 


143.75 ±0.67 


143.70± 0.44 


144.35 ± 0.53 


144.73 ± 0.62 



Table 6. NLO cross sections for pp — > bbbb + X at the LHC (\/s = 14 TeV) in comparison 



[8] 

'NLO 



, and a corrected result based on the calculation of Ref. [8], 



with published results of Ref. [8], 

^NLO (private communication) . Our results are shown for two different subtraction schemes, 

the Catani- Seymour (CS) dipole subtraction, without (cimax = 1 j o,nd with (amax = 0.01 j restriction 
on the phase space of the subtraction, and the new Nagy-Soper (NS) scheme, including the numerical 
error from the Monte Carlo integration. The scale choice is fiR = pp = fiQ, with po ^s defined in 
Eq. (3.1), and the cross sections are evaluated with the CTEQ6.5 pdf set. 



Also note that in Ref. [8] an NLO pdf set has been used both for the LO and NLO 
result. Although this may be justified to study the impact of higher-order corrections 
to the partonic cross section, it can be misleading when establishing the genuine effect 
of higher-order corrections, the adequacy of the scale choice or the size/shape of inte- 
grated/differential i^-factors used in experimental analyses when comparing Monte Carlo 
simulations to the LHC data. 
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4 Summary 

The production of four bottoms quarks, pp — )■ bbbb +X, provides an important background 
for new physics searches at the LHC. We have performed a calculation of the NLO QCD 
corrections to this process with the Helac-NLO framework, employing a new subtraction 
scheme for treating real radiation corrections at NLO implemented in Helac-Dipoles. 
Results have been presented for inclusive and differential cross-sections for pp — >• bbbb + X 
at the LHC at the centre-of-mass energy of ^/s = 14 TeV. We find that the higher-order 
corrections significantly reduce the scale dependence, with a residual theoretical uncertainty 
of about 30% at NLO. The impact of the bottom quark mass is moderate for the cross 
section normalisation and negligible for the shape of distributions. The fully differential 
NLO cross section calculation for the process pp — )• bbbb + X presented in this paper 
provides an important input for the experimental analyses and the interpretation of new 
physics searches at the LHC. 
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